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Abstract — Furnaces are most commonly used for melting 
of materials. Induction furnaces are more beneficial as no 
fuel is required. It is a problem to find life cycle of 
Induction Melting Furnace Wall under load variation. The 
induction melting furnace wall is made of alumina ramming 
mass which is one kind of refractory material. The failure 
occurs due to cyclic thermal stresses due to heating and 
cooling cycles. Temperature distribution and thermal stress 
distribution fields of the induction melting furnace 
refractory wall were calculated by using explicit finite 
difference analysis based on the physical description of its 
failure under low cycle thermal fatigue conditions. The life 
span of the refractory wall is required to be found out by 
means of critical thermal stresses created inside the 
refractory wall of induction melting furnace wall from 
modified S - log N Curve. 

Keywords — Advanced stress life analysis, Temperature 
distribution, Stress distribution, explicit finite difference 
method, Alumina ramming mass. 

I. INTRODUCTION 

Furnace is a term used to identify a closed space here heat is 
applied to a body in order to raise its temperature. The 
source of heat may be fuel or electricity. Commonly, metals 
and alloys and sometimes non-metals are heated in 
furnaces. The purpose of heating defines the temperature of 
heating and heating rate. Increase in temperature softens the 
metals. They become amenable to deformation. This 
softening occurs with or without a change in the metallic 
structure. Heating to lower temperatures (below the critical 
temperature) of the metal softens it by relieving the internal 
stresses. On the other hand, metals heated to temperatures 
above the critical temperatures leads to changes in crystal 
structures and recrystallization like annealing. Further some 
metals and alloys are melted, ceramic products vitrified, 
coals coked, metals like zinc are vaporized and many other 
processes are performed in Furnaces. 

Induction furnaces are widely used in the iron industry for 
the casting of the different grades of cast iron products. 


Refractory wall of induction melting furnace is a key 
component which is used as insulation layer. It is made of 
ramming mass like silica, alumina, magnesia etc. The 
refractory wall is directly influenced by the thermal cycling 
of the high temperature molten iron in the furnace. Thermal 
fatigue failure is easy to happen for it because of the larger 
phase transformation thermal stresses and it has a shorter 
life. This can cause serious production accidents. Therefore, 
the service life problem of the refractory wall has always 
been a focus of attention in the application of this to the 
industry. (A V K Suryanarayana, Fuels Furnaces Refractory 
and Pyrometry). 

A comprehensive literature review on computational 
investigation on different kinds of furnaces is done to study 
research trends. (Nirajkumar Mehta, May 2012). A review 
is done on applications of different numerical methods in 
heat transfer with its applications. (N C Mehta, Vipul B 
Gondaliya et al. February 2013). Thermal fatigue analysis 
of induction melting furnace wall is done for silica ramming 
mass. (N C Mehta, Akash D Raiyani et al. February 2013). 
A review is done for research on induction heating. (Vimal 
R Nakum et al, April 2013). A review is done for metal 
forming analysis using different numerical methods. (N C 
Mehta et al. May 2013). Transient heat transfer analysis of 
induction furnace is done by using finite element analysis. 
(Vipul Gondaliya et al, August 2013). Thermal fatigue 
analysis of induction furnace wall is done for alumina 
ramming mass. (N C Mehta et al, October 2013). 

Thermal analysis of hot wall condenser is done for domestic 
refrigerator using numerical method for temperature 
distribution. (Akash D Raiyani et al, July 2014). 
Optimisation of wall thickness is done for minimum heat 
loss for induction furnace by finite element analysis. 
(Dipesh D Shukla et al, December 2014). A review is done 
on numerical analysis of furnace. (N C Mehta et al, April 
2015). Thermal fatigue analysis of induction furnace wall is 
done for zirconia. (Nirajkumar C Mehta et al, April 2015). 
Comparison of finite difference method and finite element 
method is done for 2 D transient heat transfer problem. 
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(Nirajkumar Mehta et al, April 2015). Thermal fatigue 
analysis of induction furnace wall is done for magnesia 
ramming mass. (Nirajkumar C Mehta et al, June 2015). 
Advanced mathematical modelling of heat transfer is done 
for induction furnace wall of zirconia. (Nirajkumar C Mehta 
et al, December 2016). Advanced heat transfer analysis is 
done for alumina based refractory wall of induction furnace. 
(Nirajkumar Mehta et al, December 2016). 

Here, Explicit Finite Difference Method is used to find out 
temperature and thermal stress variation with respect to 
time. 

H. DEVELOPMENT OF ADVANCED EXPLICIT 
FINITE DIFFERENCE MODEL 

We have divided Induction Furnace Wall into a Nodal 
Network as shown in Fig. 1. It is divided into 24 nodes. We 
have derived Explicit Finite Difference Equations for all 
nodes as per the boundary conditions applied to it. The 
furnace wall is having thermal conduction heat transfer 
between different nodes. It is having atmospheric heat 
convection ha applied from top side of the furnace wall 
which is open to atmosphere. It is having heat convection 
from molten metal from inside which is hi. It is having heat 
convection ho from cooling water which is circulating 
outside the furnace wall. (Yunus A Cengel, Heat and Mass 
Transfer). 

To solve this advanced heat transfer problem of induction 
melting furnace wall which is made from Alumina 
Ramming Mass, the following initial and boundary 
conditions, material properties and basic assumptions are 
made: 

Refractory Materials for induction furnace wall meets the 
basic assumptions in the science of mechanics. 

> Environmental Temperature is homogeneous at 
27° C. 

r Ignore the influence of heat radiation. 

> Ignore the effect of gravity field. 

> The surface of induction melting furnace wall is 
clean. 

> The initial temperature of the induction melting 
furnace is set 27° C and it is agreement with the 
ambient temperature during solving the problem. 

> Heat convections are considered constant for this 
analysis. 

> Scarp material input inside furnace is considered 
uniform for our analysis. 



Fig.l: Nodal network for finite difference method 
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III. PROGRAMMING & SOLUTION 

With the help of a computer program we can solve the 
matrix created by finite difference equations for 24 nodes. 
We can calculate temperature distribution and stress 
distribution with respect to time. 


Table. 1: Material Property and Boundary Conditions 


Material Properties and Boundary 

Conditions for Alumina Ramming Mass 

Unit 

1 

Internal Film Co-efficient 

hi 

200 

W/m 2 

K 

2 

External Film Co-efficient 

ho 

40 

W/m 2 

K 

3 

Atmosphere Film Co- 
efficient ha 

10 

W/m 2 

K 

4 

Density 

3400 

Kg/m 3 

5 

Time Interval At 

10 

Secon 
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ds 

6 

Thermal Conductivity k 

2.6 

W/m 

K 

7 

Temperature outside 

Furnace Wall 

303 

Kelvin 

8 

Temperature inside 

Furnace Wall 

1873 

Kelvin 

9 

Temperature of Air 

303 

Kelvin 

10 

Specific Heat 

920 

J/kgK 

11 

Elasticity Constant 

220000 

N/ m 2 

12 

Thermal Expansion Co- 
efficient 

0.00000 

088 

m/K 

13 

Ultimate Stress 

500 

MPa 


Computer Program 

#include<iostream.h> 

#include<conio . h> 

#include<math.h> 

#include<complex.h> 

#include<fstream. h> 
void main() 

{ clrscr; 

ofstreammyfile; 

myfile.open ("outputfdm.txt"); 

myfile«" Writing this to a fileAn"; 

float T[25] [500]; 

floats [25] [500]; 

float hi,ho,ha,To,r,t,k,Ta,c,Th,E,Af,Sut,Se,m,N,L; 
cout«"Enter in W/m2 K value of hi"; 
cin»hi; 

cout«"Enter in W/m2 K value of ho"; 
cin»ho; 

cout«"Enter in W/m2 K value of ha"; 
cin»ha; 

cout«"Enter in Kg/m3 value of density"; 
cin»r; 

cout«"Enter in seconds time interval delta t"; 
cin»t; 

cout«"Enter in W/m K value of thermal conductivity"; 
cin»k; 

cout«"Enter in Kelvin Temperature outside Furnace 

Wall"; 

cin»Ta; 

cout«"Enter in Kelvin Temperature inside Furnace Wall"; 
cin»Th; 

cout«"Enter in Kelvin Temperature of Air"; 
cin»To; 

cout«"Enter in J/Kg K value of Specific Heat"; 
cin»c; 


cout«"Enter in N/m2 value of Elasticity Constant"; 
cin»E; 

cout«"Enter in m/K Thermal Expansion Co-efficient"; 
cin»Af; 

cout«"Enter in MPA value of Ultimate Stress"; 
cin»Sut; 

T[l][0]=300; 

T[2][0]=300; 

T[3][0]=300; 

T[4][0]=300; 

T[5][0]=300; 

T[6][0]=300; 

T[7][0]=300; 

T[8][0]=300; 

T[9][0]=300; 

T[10][0]=300; 

T[ll][0]=300; 

T[12][0]=300; 

T[13][0]=300; 

T[14][0]=300; 

T[15][0]=300; 

T[16][0]=300; 

T[17][0]=300; 

T[18][0]=300; 

T[19][0]=300; 

T[20][0]=300; 

T[21][0]=300; 

T[22][0]=300; 

T[23][0]=300; 

T[24][0]=300; 

for (inti=0; i<=270;i++) 


T[l][i+1] = (((0.5*ha*0.06*(To-T[l][i])) +(ho*0.3*(Ta- 
T[ 1 ] [i])/2)+(0.5*k*0.3*(T[2] [i]- 
T[l][i])/0.06)+(0.5*k*0.06*(T[4][i]-T[l][i])/0.3) 
)*((4*t)/(r*c*0.06*0.3)))+T[l] [i] ; 

T[2][i+1] = (((ha*0.06*(To-T[2][i]))+(0.5*k*0.3*(T[l][i]- 
T[2][i])/0.06)+(0.5*k*0.3*(T[3][i]-T[2][i])/0.06) 
+(0.5*k*0.06*(T[5][i]- 
T[2][i])/0.3))*((2*t)/(r*c*0.06*0.3)))+T[2][i]; 

T[3][i+1] = (((ha*0.06*(To-T[3][i])*0.5) +( hi*0.3*(Th- 
T[3][i])*0.5 ) +(0.5*k*0.3*(T[2][i]- 

T[3][i])/0.06)+(0.5*k*0.06*(T[6][i]-T[3][i])/0.3)) 
*((4*t)/(r*c*0.06*0.3)))+T[3][i]; 

T[4][i+1] = ((( ho*0.3*(Ta-T[4][i])) +(k*0.3*(T[5][i]- 
T[4] [i])/0.06)+(0.5*k*0.06*(T[ 1 ] [i] -T[4] [i])/0.3) + 

(0.5*k*0.06*(T[7][i]-T[4][i])/0.3)) 
*((2*t)/(r*c*0.06*0.3)))+T[4][i]; 
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T[5][i+1] = (((k*0.3*(T[4][i]- 

T[5] [i])/0.06)+(k*0.3*(T[6] [i] -T[5] [i])/0.06) + 

(k*0.06*(T[8][i]-T[5][i])/0.3)) 
*((t)/(r*c*0.06*0.3)))+T[5][i]; 

T[6][i+1] = (((hi*0.3*(Th-T[6][i])) + (k*0.3*(T[2][i]- 
T[6][i])/0.06)+(0.5*k*0.06*(T[3][i]-T[6][i])/0.3) + 

(0.5*k*0.06*(T[9][i]-T[6][i])/0.3)) 
*((2*t)/(r*c*0.06*0.3)))+T[6][i]; 

T[7][i+1] = (((ho*0.3*(Ta-T[7][i])) + (k*0.3*(T[8][i]- 
T[7][i])/0.06)+(0.5*k*0.06*(T[4][i]-T[7][i])/0.3) + 

(0.5*k*0.06*(T[10][i]-T[7][i])/0.3)) 
*((2*t)/(r*c*0.06*0.3)))+T[7][i]; 

T[8][i+1] = (((k*0.3*(T[7][i]-T[8][i])/0.06)+( 

k*0.3 *(T[9] [i]-T[8] [i])/0.06 ) +(k*0.06*(T[5][i]- 

T[8][i])/0.3) + (k*0.06*(T[ 1 1 ] [i] -T[8][i])/0.3)) 

*((t)/(r*c*0.06*0.3)))+T[8][i]; 

T[9][i+1] = (( (hi*0.3*(Th-T[9][i]))+(k*0.3*(T[8][i]- 

T[9][i])/0.06)+ (0.5*k*0.06*(T[6][i]-T[9][i])/0.3) + 

(0.5*k*0.06*(T[12][i]-T[9][i])/0.3)) 
*((2*t)/(r*c*0.06*0.3)))+T[9][i]; 

T[10][i+1] = (( (ho*0. 1 8*(Ta-T[10] [i]))+(k*0. 1 8 *(T[ 11] [i] - 
T[10][i])/0.06)+ (0.5*k*0.06*(T[7][i]-T[10][i])/0.18) + 
(0.5*k*0.06*(T[15][i]-T[10][i])/0.18)) 

*((2*t)/(r*c*0.06*0. 1 8)))+T[ 10] [i] ; 

T[1 l][i+l] = (( (k*0.18*(T[10][i]-T[ll][i])/0.06)+ 

(k*0.18*(T[12][i]-T[ll][i])/0.06)+ (k*0.06*(T[8][i]- 

T[ 11] [i])/0. 1 8) + (k*0.06*(T[ 16] [i] -T[ 1 1 ] [i])/0. 18)) 

*((t)/(r*c *0.06 *0. 1 8)))+T[ 1 1 ] [i] ; 

T[12][i+1] = (( (0.5*hi*0.12*(Th-T[12][i]))+( 

0.5 *hi*0. 1 8 *(Th-T[ 12] [i])) +(k*0. 1 8*(T[ 1 1 ] [i] - 

T[12][i])/0.12)+ (0.5*k*0.18*(T[13][i]- 

T[ 1 2] [i])/0. 12)+(0.5 *k*0. 12*(T[9] [i] -T[12] [i])/0. 1 8) + 

(k*0.12*(T[17][i]-T[12][i])/0.18)) 

*((4*t)/(3 *r*c*0. 1 2*0. 1 8))) +T[ 1 2] [i] ; 

T[13][i+1] = (( (hi*0.18*(Th- 

T[13][i]))+(0.5*k*0.18*(T[12][i]-T[13][i])/0.18)+ 
(0.5*k*0.18*(T[14][i]-T[13][i])/0.18)+ (k*0.18*(T[18][i]- 
T[ 1 3] [i])/0. 18)) *((2*t)/(r*c*0. 1 8*0. 1 8)))+T[ 1 3] [i]; 
T[14][i+1] = (( (0.5*hi*0.18*(Th- 

T[ 14] [i]))+(0.5*hi*0. 1 8*(Th-T[ 14] [i])) 
+(0.5*k*0.18*(T[13][i]-T[14][i])/0.18)+ 
(0.5*k*0.18*(T[19][i]-T[14][i])/0.18)) 

*((4*t)/(r*c *0. 1 8 *0. 1 8)))+T[ 14] [i] ; 

T[15][i+1] = (( (ho*0.06*(Ta-T[15][i]))+(k*0.06*(T[16][i]- 
T[15][i])/0.06)+ (0.5*k*0.06*(T[10][i]-T[15][i])/0.06)+( 

0.5*k*0.06*(T[20][i]-T[15][i])/0.06)) 
*((2*t)/(r*c*0.06*0.06)))+T[ 15] [i] ; 


[Vol-3, Issue-1, Jan-Feb, 2017] 
ISSN: 2455-5304 

T[16][i+1] = (( (k*0.06*(T[15][i]-T[16][i])/0.06)+ 

(k*0.06*(T[ 17] [i]-T[ 16] [i])/0.06 )+(k*0.06*(T[ 1 1 ] [i]- 

T[ 16] [i])/0.06)+( k*0.06*(T[2 1 ] [i] -T[ 16] [i])/0.06)) 

*((t)/(r *c*0.06*0.06)))+T[ 1 6] [i] ; 

T[17][i+1] = (( (k*0.06*(T[ 1 6] [i] -T[ 1 7][i])/0.12)+ 

(k*0.06*(T[ 1 8] [i] -T[ 17] [i])/0. 12 )+(k*0.12*(T[12][i]- 

T[17][i])/0.06)+( k*0.12*(T[22][i]-T[17][i])/0.06)) 

*((t)/(r*c*0.12*0.06)))+T[17][i]; 

T[18][i+1] = (( (k*0.06*(T[ 17] [i] -T[ 1 8] [i])/0. 18)+ 

(k*0.06*(T[ 19] [i]-T[18] [i])/0. 18 )+(k*0.18*(T[13][i]- 

T[ 1 8] [i])/0.06)+f k*0. 1 8*(T[23] [i] -T[ 18] [i])/0.06)) 

*((t)/(r*c*0. 1 8*0.06)))+T[ 1 8] [i] ; 

T[19][i+1] = (( (hi*0.06*(Th-T[ 19] [i]))+(k*0.06*(T[ 1 8] [i]- 
T[19][i])/0.18)+ (0.5*k*0.18*(T[14][i]-T[19][i])/0.06)+( 

0.5 *k*0. 1 8 *(T[24] [i] -T[ 1 9] [i])/0.06)) 

*((2*t)/(r*c *0. 1 8 *0.06)))+T[ 1 9] [i] ; 

T[20][i+1] = (( (0.5*ho*0.06*(Ta-T[20][i]))+( 

0.5*ho*0.06*(Ta-T[20][i]))+(0.5*k*0.06*(T[21][i]- 
T[20][i])/0.06)+ (0.5*k*0.06*(T[ 15] [i]- 

T[20][i])/0.06))*((4*t)/(r*c*0.06*0.06)))+T[20][i]; 

T[21][i+1] = (( (0.5*ho*0.06*(Ta- 

T[2 1 ] [i]))+(0.5*k*0.06*(T[20] [i] -T[2 1 ] [i])/0.06)+( 
0.5*k*0.06*(T[22] [i] -T[2 1] [i])/0.06)+ (k*0.06*(T[ 16] [i]- 

T[2 1 ] [i])/0.06))*((2*t)/(r*c*0.06*0.06)))+T[2 1 ] [i] ; 
T[22][i+1] = (( (ho*0.12*(Ta- 

T[22] [i]))+(0.5*k*0.06*(T[2 1 ] [i] -T[22] [i])/0. 1 2)+( 
0.5*k*0.06*(T[23][i]-T[22][i])/0.12)+ (k*0. 12*(T[ 17] [i]- 

T[22][i])/0.06))*((2*t)/(r*c*0.12*0.06)))+T[22][i]; 
T[23][i+1] = (( (ho*0.18*(Ta- 

T[23][i]))+(0.5*k*0.06*(T[22][i]-T[23][i])/0.18)+( 

0.5 *k*0.06*(T[24] [i] -T[23] [i])/0. 1 8 )+ (k*0. 1 8 *(T[ 1 8] [i] - 

T[23][i])/0.06))*((2*t)/(r*c*0. 1 8*0.06)))+T[23] [i] ; 
T[24][i+1] = (( (0.5*ho*0. 1 8*(Ta-T[24] [i]))+( 

0.5*hi*0.06*(Th-T[24][i])) +(0.5*k*0.06*(T[23][i]- 

T[24] [i] )/0. 18)+ (0.5 *k*0. 1 8 *(T[ 1 9] [i] - 

T[24] [i])/0.06))*((4*t)/(r*c*0. 1 8*0.06)))+T[24] [i] ; 

S [ 14] [i] = E*Af*T[14][i]; 

} 

for (i=271; i<=360;i++) 

{ 

hi=ha; 

Th=Ta; 

T[l][i+1] = (((0.5*ha*0.06*(To-T[l][i])) +(ho*0.3*(Ta- 
T[ 1 ] [i])/2)+(0.5*k*0.3*(T[2] [i] - 
T[ 1 ] [i])/0.06)+(0.5*k*0.06*(T[4] [i] -T[ 1 ] [i])/0.3) 
)*((4*t)/(r*c*0.06*0.3)))+T[l][i]; 

T[2][i+1] = (((ha*0.06*(To-T[2][i]))+(0.5*k*0.3*(T[l][i]- 
T[2][i])/0.06)+(0.5*k*0.3*(T[3][i]-T[2][i])/0.06) 
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+(0.5*k*0.06*(T[5][i]- 

T[2][i])/0.3))*((2*t)/(r*c*0.06*0.3)))+T[2][i]; 

T[3][i+1] = (((ha*0.06*(To-T[3][i])*0.5) +( hi*0.3*(Th- 
T[3][i])*0.5 ) +(0.5*k*0.3*(T[2][i]- 

T[3][i])/0.06)+(0.5*k*0.06*(T[6][i]-T[3][i])/0.3)) 
*((4*t)/(r*c*0.06*0.3)))+T[3][i]; 

T[4][i+1] = ((( ho*0.3*(Ta-T[4][i])) +(k*0.3*(T[5][i]- 
T[4][i])/0.06)+(0.5*k*0.06*(T[l][i]-T[4][i])/0.3) + 

(0.5*k*0.06*(T[7][i]-T[4][i])/0.3)) 
*((2*t)/(r*c*0.06*0.3)))+T[4][i]; 

T[5][i+1] = (((k*0.3*(T[4][i]- 

T[5] [i])/0.06)+(k*0.3*(T[6] [i] -T[5] [i])/0.06) + 

(k*0.06*(T[8][i]-T[5][i])/0.3)) 
*((t)/(r*c*0.06*0.3)))+T[5][i]; 

T[6][i+1] = (((hi*0.3 *(Th-T[6] [i])) + (k*0.3*(T[2][i]- 
T[6][i])/0.06)+(0.5*k*0.06*(T[3][i]-T[6][i])/0.3) + 

(0.5*k*0.06*(T[9] [i] -T[6] [i])/0.3)) 
*((2*t)/(r*c*0.06*0.3)))+T[6][i]; 

T[7][i+1] = (((ho*0.3*(Ta-T[7][i])) + (k*0.3*(T[8][i]- 
T[7][i])/0.06)+(0.5*k*0.06*(T[4][i]-T[7][i])/0.3) + 

(0.5 *k*0.06*(T[ 1 0] [i]-T[7] [i])/0.3)) 
*((2*t)/(r*c*0.06*0.3)))+T[7][i]; 

T[8][i+1] = (((k*0.3 *(T[7] [i] -T[8] [i])/0.06)+( 

k*0.3*(T[9][i]-T[8][i])/0.06 ) +(k*0.06*(T[5][i]- 

T[8][i])/0.3) + (k*0.06*(T[ll][i]-T[8][i])/0.3)) 

*((t)/(r*c*0.06*0.3)))+T[8][i]; 

T[9][i+1] = (( (hi*0.3*(Th-T[9][i]))+(k*0.3*(T[8][i]- 

T[9][i])/0.06)+ (0.5*k*0.06*(T[6][i]-T[9][i])/0.3) + 

(0.5*k*0.06*(T[12][i]-T[9][i])/0.3)) 
*((2*t)/(r*c*0.06*0.3)))+T[9][i]; 

T[10][i+1] = (( (ho*0. 1 8*(Ta-T[ 1 0] [i]))+(k*0. 1 8 *(T[ 1 1 ] [i] - 
T[10][i])/0.06)+ (0.5 *k*0.06*(T[7] [i]-T[ 1 0] [i])/0. 1 8) + 

(0.5*k*0.06*(T[15][i]-T[10][i])/0.18)) 

*((2*t)/(r*c*0.06*0. 1 8)))+T[ 10] [i] ; 

T[ll][i+1] = (( (k*0.18*(T[10][i]-T[ll][i])/0.06)+ 

(k*0.18*(T[12][i]-T[ll][i])/0.06)+ (k*0.06*(T[8][i]- 

T[1 l][i])/0. 18) + (k*0.06*(T[ 1 6] [i] -T[ 1 1 ] [i] )/0. 1 8)) 

*((t)/(r*c*0.06 *0. 1 8)))+T[ 1 1 ] [i] ; 

T[12] [i+1] = (( (0.5*M*0.12*(Th-T[12][i]))+( 

0.5 *hi*0. 1 8*(Th-T[ 12] [i])) +(k*0.18*(T[ll][i]- 

T[12][i])/0.12)+ (0.5*k*0.18*(T[13][i]- 

T[ 1 2] [i])/0. 12)+(0.5 *k*0. 12*(T[9] [i] -T[12] [i])/0. 1 8) + 

(k*0.12*(T[17][i]-T[12][i])/0.18)) 

*((4*t)/(3 *r*c *0. 1 2*0. 1 8))) +T[ 1 2] [i] ; 

T[13][i+1] = (( (hi*0.18*(Th- 

T[ 1 3] [i]))+(0.5 *k*0. 1 8*(T[ 12] [i] -T[ 1 3] [i])/0. 18)+ 
(0.5*k*0.18*(T[14][i]-T[13][i])/0.18)+ (k*0.18*(T[18][i]- 
T[ 1 3] [i])/0. 18)) *((2*t)/(r*c*0. 1 8*0. 1 8)))+T[ 1 3] [i] ; 
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T[ 14] [i+1] = ~ (0.5*hi*0.18*(Th- 

T[14][i]))+(0.5*hi*0.18*(Th-T[14][i])) 

+(0.5*k*0.18*(T[13][i]-T[14][i])/0.18)+ 

(0.5*k*0.18*(T[19][i]-T[14][i])/0.18)) 

*((4*t)/(r*c *0. 1 8 *0. 1 8)))+T[ 14] [i] ; 

T[15][i+1] = (( (ho*0.06*(Ta-T[15][i]))+(k*0.06*(T[16][i]- 
T[15][i])/0.06)+ (0.5 *k*0.06*(T[ 10] [i]-T[ 1 5] [i])/0.06)+( 

0.5*k*0.06*(T[20] [i] -T[ 15] [i])/0.06)) 
*((2*t)/(r*c*0.06*0.06)))+T[15][i]; 

T[16][i+1] = (( (k*0.06*(T[15][i]-T[16][i])/0.06)+ 

(k*0.06*(T[ 17] [i] -T[ 16] [i])/0.06 )+(k*0.06*(T[ 1 1 ] [i]- 

T[ 16] [i])/0.06)+( k*0.06*(T[2 1 ] [i] -T[ 16] [i])/0.06)) 

*((t)/(r*c*0.06*0.06)))+T[16][i]; 

T[17][i+1] = (( (k*0.06*(T[16][i]-T[17][i])/0.12)+ 

(k*0.06*(T[18][i]-T[17][i])/0.12 )+(k*0.12*(T[12][i]- 

T[ 17] [i])/0.06)+( k*0.12*(T[22] [i] -T[ 17] [i])/0.06)) 

*((t)/(r*c*0.12*0.06)))+T[17][i]; 

T[18][i+1] = (( (k*0.06*(T[17][i]-T[18][i])/0.18)+ 

(k*0.06*(T[ 19] [i]-T[ 1 8] [i])/0. 1 8 )+(k*0. 1 8*(T[ 1 3] [i]- 

T[ 1 8] [i])/0.06)+( k*0. 1 8*(T[23] [i] -T[ 18] [i])/0.06)) 

*((t)/(r*c*0. 1 8*0.06)))+T[ 1 8] [i] ; 

T[19][i+1] = (( (hi*0.06*(Th-T[ 19] [i]))+(k*0.06*(T[ 1 8] [i]- 
T[19][i])/0.18)+ (0.5*k*0.18*(T[14][i]-T[19][i])/0.06)+( 

0.5 *k*0. 1 8 *(T[24] [i] -T[ 1 9] [i])/0.06)) 

*((2*t)/(r*c *0.18 *0.06)))+T[ 1 9] [i] ; 

T[20][i+1] = (( (0.5*ho*0.06*(Ta-T[20][i]))+( 

0.5*ho*0.06*(Ta-T[20][i]))+(0.5*k*0.06*(T[21][i]- 
T[20] [i] )/0.06)+ (0.5 *k*0.06*(T[ 1 5] [i] - 

T[20][i])/0.06))*((4*t)/(r*c*0.06*0.06)))+T[20][i]; 
T[21][i+1] = (( (0.5*ho*0.06*(Ta- 

T[2 1 ] [i]))+(0.5*k*0.06*(T[20] [i] -T[2 1 ] [i])/0.06)+( 
0.5*k*0.06*(T[22] [i] -T[2 1] [i])/0.06)+ (k*0.06*(T[ 16] [i]- 

T[21] [i])/0.06))*((2*t)/(r*c*0.06*0.06)))+T[2 1 ] [i] ; 
T[22][i+1] = (( (ho*0.12*(Ta- 

T[22] [i]))+(0.5*k*0.06*(T[2 1 ] [i] -T[22] [i])/0. 1 2)+( 

0.5 *k*0.06 *(T[23] [i] -T[22] [i])/0. 1 2 )+ (k*0. 12*(T[ 17] [i]- 

T[22][i])/0.06))*((2*t)/(r*c*0.12*0.06)))+T[22][i]; 
T[23][i+1] = (( (ho*0.18*(Ta- 

T[23] [i]))+(0.5*k*0.06*(T[22] [i]-T[23] [i])/0. 1 8)+( 

0.5 *k*0.06*(T[24] [i] -T[23] [i])/0. 1 8 )+ (k*0. 1 8 *(T[ 1 8] [i] - 

T[23] [i])/0.06))*((2*t)/(r*c*0. 1 8*0.06)))+T[23] [i] ; 
T[24][i+1] = (( (0.5*ho*0. 18*(Ta-T[24] [i]))+( 

0.5*hi*0.06*(Th-T[24][i])) +(0.5*k*0.06*(T[23][i]- 

T[24] [i] )/0. 18)+ (0.5 *k*0. 1 8*(T[ 1 9] [i] - 

T[24] [i] )/0.06))*((4*t)/(r*c *0. 1 8 *0.06)))+T[24] [i] ; 

S [ 14] [i] = E*Af*T[14][i]; 

} 

for(i=0; i<=360;i=i+6) 
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myfile«"MAXIMUM TEMPERATURE AFTER 
"«i/6«" MINUTES"«T[14][i]«" KELVIN"«endl; 
myfile«"MAXIMUM STRESS AFTER "«i/6«" 

MINUTES "«S[14][i]«" MPa"«endl; } 

for(i=0; i<=300; i++) 

{ 

if (i==270) 

{ 

myfile«" S "«S [14] [i] ; 

Se=0.15*Sut; 
m = (Sut-Se)/7; 

N = 7 -((S[14][i]-Se)/m); 

L = pow(10,N); 
myfile«"N"«N; 
myfile«"LIFE CYCLE"«L; 
myfile.close(); 

} 

} 

getch(); } 

IV. RESULTS AND DISCUSSION 

We can see from the Fig. 2 that maximum temperature is 
increasing from atmospheric temperature 300 K and reaches 
to maximum temperature 1787 K in 45 minutes and then 
starts reducing and reaches to 869 K in next 15 minutes. It 
again starts increasing and reaches to maximum 1787 K 
after 105 minutes and again starts reducing. There are 10 
similar temperature cycles in one day. 


Maximum Temperature v/s Time 
Alumina Ramming Mass 

2000 



0 

Time (Minutes) 

Fig. 2: Maximum temperature v/s time graph for alumina 
ramming mass 
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Maximum Stress v/s Time Alumina 
Ramming Mass 
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Fig. 3: Maximum thermal stress v/s time graph for alumina 
ramming mass 

We can see from the Fig. 3 that maximum thermal stress is 
increasing from initial condition 0 MPa and reaches to 
maximum stress 346 MPa in 45 minutes and then starts 
reducing and reaches to 168 MPa in next 15 minutes. It 
again starts increasing and reaches to maximum stress 346 
MPa after 105 minutes and again it starts reducing. There 
are 10 similar thermal stress cycles in one day. 

Stress-Life Method 

To determine life of any component by Stress-Life Method, 
we need to find out ultimate strength and endurance limit of 
the component for the required material. 

We know the values of ultimate stress for these all 
materials. Alumina ramming mass is having ultimate 
strength of 500 MPa. We can find out Se’ from the 
equation given below or we can say dividing value of 
ultimate strength. 

We know value of ultimate stress of alumina 
ramming mass is 500 Mpa. 

We know the relation between Sut and Se’ so that we can 
find out Se’. 

Se’ = 0.5 * Sut= 250 MPa 
Endurance Limit Modifying Factors 

We have seen that the rotating-beam specimen used in the 
laboratory to determine endurance limits is prepared very 
carefully and tested under closely controlled conditions. 

It is unrealistic to expect the endurance limit of a 
mechanical or structural member to match the values 
obtained in the laboratory. 
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Stress v/s Life Cycle Alumina Ramming 
Mass 


Some differences include 

Material: composition, basis of failure, variability 
Manufacturing: method, heat treatment, fretting corrosion, 
surface condition, stress concentration 
Environment: corrosion, temperature, stress state, relaxation 
times 

Design: size, shape, life, stress state, stress concentration, 
speed, fretting, galling 

Marin identified factors that quantified the effects of surface 
condition, size, loading, temperature, and miscellaneous 
items. The question of whether to adjust the endurance limit 
by subtractive corrections or multiplicative corrections was 
resolved by an extensive statistical analysis of a 4340 
(electric furnace, aircraft quality) steel, in which a 
correlation coefficient of 0.85 was found for the 
multiplicative form and 0.40 for the additive form. 

A Marin equation is therefore written as 

Se = ka*kb*kc*kd*ke*kf* Se’ 

Where, 

ka= surface condition modification factor = 0.84 
kb= size modification factor = 0.98 
kc= load modification factor = 0.96 
kcl= temperature modification factor = 0.45 
ke= reliability factor = 0.9 

kf= miscellaneous-effects modification factor = 0.95 
Se’= specimen endurance limit 

Se = endurance limit at the critical location of a machine 
part in the geometry and condition of use. 

We can find out different factor like surface finish factor, 
size factor, loading factor, temperature factor, reliability 
factor, miscellaneous effects factor as per the guideline. 
(Joseph E. Shigleyef al. Machine Engineering Design) 

Now, we can find out endurance limit for all different 
materials. 

Se = ka*kb*kc*kd*ke*kf*Se' = 75 MPa 


600 
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Life Cycle (Log N) 

Fig. 3: Stress v/s Life Cycle for Alumina Ramming Mass 

N (No. of working cycles) = L n ii n X (100-Probability) / (100 
X Miscellaneous Factor) 

D (No. of working days) = N / (No. of working cycles per 
day) 


Table. 2: Number of Working Cycles and Number of 
Working Days with Probability 


Sr. 

No. 

Probability 

(%) 

No. of working 
cycles for 

Alumina 

Ramming Mass 

No. of 

days 

working 

1 

100 

342 

38.0 

2 

90 

349 

38.8 

3 

80 

356 

39.5 

4 

70 

363 

40.3 

5 

60 

369 

41.0 

6 

50 

376 

41.8 

7 

40 

383 

42.6 

8 

30 

390 

43.3 

9 

20 

397 

44.1 

10 

10 

404 

44.8 


We can understand from this table that all cases of 
induction furnace wall will work for minimum 342 working 
cycles and 38 days. The probability will be decreasing 
continuously with increase in life span and working days. It 
happens very rarely that induction furnace wall will work 
for 404 working cycles and 44.8 days so its probability is 
minimum 10% of cases can be having such a large life span. 
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V. CONCLUSION 

Induction Melting Furnaces are highly used now- a-days for 
melting of different kinds of materials. The problem comes 
from the alumina ramming mass of losing its material 
properties and failure occurs within 340-370 hours of 
lifetime. It will disturb production schedule as it requires 
time to replace the induction melting furnace wall of 
alumina ramming mass. Explicit finite difference analysis is 
done for induction melting furnace refractory wall and 
validation is done with respect to experimental Results. 
Explicit finite difference analysis is done with respect to 
actual working conditions of induction furnace and material 
properties of alumina ramming mass. Then S - log N 
Curves are plotted for Life Span Prediction. 

We had found life span of different materials like alumina 
ramming mass is 340-370 cycles. From the results of 
experimental study and explicit finite difference analysis of 
thermal fatigue failure of induction melting furnace wall, it 
can be seen that finite difference model exactly predicts the 
failure of the induction furnace refractory wall and the 
definite solution conditions in the finite difference 
numerical calculation are accurate. The fatigue life of the 
induction melting furnace refractory wall under thermal 
fatigue working conditions was predicted using stress-life 
method by plotting S - log N curves for alumina ramming 
mass on the basis of explicit finite difference calculations 
and maximum thermal stress in the induction melting 
furnace refractory wall. We can use it as a linear to increase 
lifespan or we can use premixed alumina ramming mass for 
economical and better working lifespan of induction 
melting furnace wall. 

The accuracy of the fatigue life prediction for the induction 
melting furnace refractory wall depends upon temperature 
and thermal stress spectrum calculated at the critical point 
by explicit finite difference method and S-log N curves 
prepared from the material properties and boundary 
condition for alumina ramming mass. 
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